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Abstract 

Using the method of anharmonic lattice statics, we calculate the equilibrium structure of steps on 
180° ferroelectric domain walls (DW) in PbTiOa. We consider three different types of steps: i) Ti- 
Ti step that joins a Ti-centered DW to a Ti-centered DW, (ii) Pb-Pb step that joins a Pb-centered 
DW to a Pb-centered DW, and (iii) Pb-Ti step that joins a Pb-centered DW to a Ti-centered DW. 
We show that atomic distortions due to these steps broaden a DW but are localized, i.e., they are 
confined to regions with dimensions of a few lattice spacings. We see that a step locally thickens 
the domain wall; the defective domain wall is two to three times thicker than the perfect domain 
wall depending on the step type. We also observe that steps distort the polarization distribution 
in a mixed Bloch-Neel like way; polarization rotates out of the domain wall plane near the steps. 
Our calculations show that Pb-Pb steps have the lowest static energy. 

PACS numbers: 75.60. Ch,77.80.Dj 



Electronic address: arash.yavari@ce.gatech.edu 



I. INTRODUCTION 



Ferroelectric materials are an important subclass of polar materials due to their wide 
range of applications in ultrasound imaging, microelectromechanical systems, high strain 
actuators, electro-optical systems, photothermal imaging, and high density storage devices 
It is known that some important properties of ferroelectric materials are due to the 
presence of domain walls, which are two-dimensional defects that separate regions with 
uniform polarization [2!] . This explains the importance of a detailed study of the properties 
of the domain walls. 

From both experimental and theoretical studies, it is observed that the thickness of 

nn nri 

domain walls can vary from a few angstroms [3H7[ to a few micrometers |8|, |9| . It has been 
suggested that this wide scatter in the domain wall thickness is due to the presence of 



point defects 
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12|. Another important property of domain walls is the behavior of the 



polarization profile near the domain wall. It is well known that 180° domain walls have an 
Ising-like nature. Using Monte-Carlo simulation, Padilla et al. \^ showed the predominant 
Ising-like character of 180° domain walls in tetragonal BaTiOs along the tetragonal axis. In 
180° domain walls, polarization vector can either rotate in a plane parallel to the domain wall 
(Bloch type) or normal to the domain wall (Neel type) [13!]. Subsequent works on the domain 
walls showed that domain walls can have mixed characters. Using density functional theory, 
Lee et al. |l^ showed that while 180° domain walls in PbTiOs are predominantly Ising-like, 
they have some Neel characters as well. Having the domain walls parallel to the (lOO)-plane, 
we know that polarization is mainly along the (OlO)-direction (see Fig. [T]). As Lee et al. [isl 
showed close to the domain wall polarization has normal components (normal to the domain 
wall) with magnitudes in the order of 1-2 percent of the bulk polarization. Angoshtari and 
Yavari 6| observed a similar behavior at finite temperatures for perfect 180° domain walls. 
They saw normal components in the order of 2 percent of the bulk polarization in their 
finite-temperature structure calculations. Recently, first-principle-based simulations have 
led to the prediction of vortex type polarization distribution in zero- dimensional ferroelectric 
nanodots Q, 1^. 

It is believed that steps have an important role in domain wall motion. Nettleton [l^ 
proposed a model for sidewise displacement of a 180° domain wall in a single crystal barium 
titanate and suggested that the formation of an irregular pattern of steps of varying shapes 
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and sizes results in the motion of the domain wall and the speed of the domain wall motion 



is determined by the rate of formation and disappearance of these steps. Shur et al. [17 1 
considered steps on 180° domain walls and proposed a mechanism for domain wall motion 
in weak and strong fields. Shin et al. pLE] used atomistic molecular dynamics and coarse- 
grained Monte Carlo simulations to analyze the nucleation and growth mechanism of domain 
walls in PbTiOa and BaTiOs. 

In this work we investigate the effect of steps, which are one- dimensional defects, on 
180° domain walls parallel to (lOO)-planes in PbTiOa using the anharmonic lattice statics 
method. We consider Ti-Ti steps that join a Ti-centered DW to another Ti-centered DW, 
Pb-Pb steps that join a Pb-centered DW to another Pb-centered DW, and Pb-Ti steps that 
join a Pb-centered DW to a Ti-centered DW. As the initial configuration, we start from the 
atomic configuration of perfect 180° domain walls and then relax the structure iteratively 
to obtain the optimized atomic configuration. 

This paper is organized as follows. In ^JTll we explain the initial geometry of steps that 
we analyzed throughout this work. In §1111 we discuss the method of anharmonic lattice 
statics and the shell potential for PbTiOs that we used in our calculations. We present our 
numerical results in §IV[ The paper ends with concluding remarks in ^|Vl 



II. GEOMETRY OF STEPS 



The geometry of the relaxed unit cell of tetragonal PbTiOs is shown in FigJTJ^a). The 
nonzero relative displacements in the (OlO)-direction between the center of the positive 
and negative charges generate a polarization in the (OlO)-direction (we are using a shell 
potential). In the 180° domain walls, direction of polarization switches across the domain 
wall. There are two types of 180° domain wall in PbTiOs, namely, Ti-centered and Pb- 
centered domain walls. Using the relaxed bulk configurations, it is possible to calculate the 
atomic structure of both types . 

Consider a domain wall parallel to a (lOO)-plane. By a step on the domain wall we mean 
the region where the domain wall joins another domain wall parallel to the first wall with 
an offset in the (lOO)-direction (see FigsiU^b) and (c)).[25] We consider three different steps: 
Ti-Ti, Pb-Pb, and Pb-Ti. FigJ2] shows the unrelaxed initial configuration for each step. Note 
that assuming that the step is limited to a single unit cell, i.e. if the two domain walls are 



3 



(c) 



DW 



DW 




e;100> I I 

I i o I 

I ►<ooi> ^ 



o 



() 



o 



step 



Case I 



Case II 



FIG. 1: (a) The relaxed configuration of the unit cell of tetragonal PbTiOs. a and c are the 
tetragonal lattice parameters. Note that 01, 02, and 03 refer to oxygen atoms located on (001), 
(100), and (OlO)-planes, respectively. 6 denotes the y-displacements of the atoms from the cen- 
terosymmetric position and arrows near each atom denote the direction of these displacements, (b) 
Schematic profile of polarization close to a step, (c) Two different possibilities for a Pb-Pb step. 



one or half a lattice spacing apart, there would be more than one possibility for the step 
configuration. As an example, we plot two possibilities for Pb-Pb step in Fig. [D^c). In this 
figure. Case I shows a Pb-Pb step in (001) PbO-plane while Case II shows another Pb-Pb 
step in (001) Ti02-plane. Note that there are still other possibilities for Pb-Pb steps. We 
should emphasize that the configurations shown in Fig. [T]^c) and Figj2]are only the initial 
configurations that we use as the starting point for finding the final equilibrium configuration. 
We observe that as far as we confine the step to a single unit cell, the anharmonic lattice 
statics iterations converge to the same solution regardless of the initial configuration of the 
step. Therefore, the exact choice of the initial step configuration is not important in the 
final equilibrium structure. We should also emphasize that we are analyzing a single step 
on a single domain wall in an infinite crystal, i.e. no periodicity assumptions are made. 
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Note that in Figj2l domain walls away from the step have polarization only along the (010)- 
direction. Note also that we assume a 2-D symmetry reduction, which means that all the 
atoms with the same x and z coordinates (x, y, and z are coordinates along the (100), (010), 
and (OOl)-directions, respectively) have the same displacements. Therefore, we partition the 
3-D lattice £ as £ = where Ciaj3 and Z are 1-D equivalence classes parallel 



to the (OlO)-direction and the set of integers, respectively. See 19|, |20| for more details on 
the symmetry reduction. 

III. METHOD OF CALCULATION 

We use the method of anharmonic lattice statics {lol to calculate the atomic structure 
of steps. We use a shell potential for PbTiOs [21] to model the atomic interactions. In this 
potential, each ion is represented by a core and a massless shell. Let C denote the collection 
of cores and shells, i C denote a core or a shell in £, and {x*},-^^ represent the current 
position of cores and shells. Then, the total static energy can be written as 



+ ^^core-shell \{^'\iizc) ' i^) 



^^short ({^^liG^) denotes short range interactions, which are assumed to be only between Pb- 
O, Ti-0, and 0-0 shells. The short range interactions are described by the Rydberg potential 
of the form {A + Br) exp(— r/C), where A, B and C are potential parameters and r is the 
distance between interacting elements. £\ong ({^'}ie£) denotes the Coulombic interactions 
between the core and shell of each ion with the cores and shells of all of the other ions. 
Note that for calculating the classical Coulombic potential and force, we use the damped 
Wolf method 3]. Finally, £'core-sheii ({x^jjg^) represents the interaction of core and shell of 
an atom and is assumed to be an anharmonic spring of the form {l/2)k2r'^ + (l/24)A;4r^ 
where ^2 and are constants. All calculations are done for absolute zero temperature. 
As is shown in Fig{T], at this temperature PbTiOa has a tetragonal unit cell with lattice 
parameters a = 3.843 A and c = 1.08a ||21]. For more details on this notation see 20]. 

For the relaxed configuration E = {x*}^^^ C R^, static energy attains a local minimum 
and hence we have 

d£ 

— = oyieC. (2) 



To obtain the solution of the above optimization problem we use the Newton method, which 
is based on a quadratic approximation near the current configuration B^: 



+ ]p^y ■Yl{B^) ■~6^ + o(\~6 



(3) 



where 6'^ = B'^^^ — B^ and H is the Hessian matrix. In the Newton method: 

5^ = -H-i {B'') ■ VS {B'') . (4) 

Having d^, the next configuration is calculated as: B''^^ = B^ + 5^. 

As the size of the simulation box increases, the calculation of the Hessian becomes inef- 
ficient and hence we use the quasi-Newton method. In this method, instead of calculating 
the Hessian in each iteration, one uses the Broyden-Fletcher-Goldfarb-Shanno (BFGS) al- 



gorithm to approximate the inverse of the Hessian 2J]. One starts with a positive-definite 



matrix and uses the following BFGS algorithm to update the Hessian at each iteration: 

C^+i = + ® - (C^ • A) ® (C- ■ A) 
{6^y - A AT . C* ■ A 

+ (A"^-C'- A)u®u, (5) 
where C' = (H')"\ A = V£'+^ - V£\ and 

C'-A 

u = — r^F — — :r- (6) 

{S'^y . A AT . C* ■ A ^ ^ 

Calculating C*"*"^, one then should use C*"*"^ instead of H^^ to update the current configura- 
tion for the next configuration B^~^^ = B^ + 6^. If C*'''^ is a poor approximation, then one 
may need to perform a linear search to refine before starting the next iteration. 

As the initial configuration for each step we start with two half lattices with the proper 
offset in the x-direction. The atomic configuration of each half lattice is the same as the 
atomic configuration in a perfect 180° domain wall. To remove the rigid body translation 
of the lattice, we fix the core of an atom in our computational box and fully relax the other 
atoms. Hence, we have SOW x (L — 3) variables in our calculations, where W and L are 
specified in FiglJl To consider the effect of the atoms outside the computational box, we 
impose rigid body translations to these atoms as the boundary conditions, i.e., we rigidly 



move the atoms outside the computational box such that they keep the perfect 180° domain 
wall configuration. To this end, we rigidly move all of the atoms outside the computational 
box in the positive (negative) direction of the z-axis with the displacements equal to the 
displacement of the first (last) atom of the first (last) row of the representative layer of 
atoms. This is marked with the red (blue) circle in Figl2]^a). The displacements of the 
atoms outside the computational box in the positive (negative) direction of the x-axis are 
equal to the displacements of the atoms in the first (last) row of the computational box that 
is located on the same column. We choose W = 20 and L = 30 as we see larger values will 
not affect the results. In all our calculations we assume force tolerance of 0.05 eV/A and 
observe that our solutions converge slowly after about 800 to 1000 iterations depending on 
the step type. Our calculations also show that using a smaller force tolerance of 0.005 eV/A 
would change the results by less than 0.1%. This justifies the above choice of the force 
tolerance. 

IV. NUMERICAL RESULTS 

In this section we present our numerical results for the three different steps as follows. As 
we mentioned earlier, x, y, and z are coordinates along the (100), (010), and (OOl)-directions, 
respectively, and the origin of the coordinate system for each step is specified in FigJU 

Pb-Pb step: The atomic configuration of Pb-cores in a Pb-Pb step are shown in FigJSt^a) 
and (b). For a clearer presentation of the atomic configuration, we have plotted the y- 
coordinates of Pb-cores for different sections v-v and h-h (see Fig|2]). In Figj3l s and d 
denote the distances of the h-h and v-v sections from the reference planes, respectively. 
Note that as is shown in Figj2l the reference plane for v-v sections (r^,) is parallel to the 
(OOl)-plane and the reference plane for h-h sections (r/j) is parallel to the (lOO)-plane. As it 
can be seen, atomic distortions in the Pb-Pb step are localized, i.e, they are confined to a 
8a X 20a box in the (OlO)-plane. Atomic distortions in the (OOl)-direction are less localized 
compared to those in the (lOO)-direction. We observe that the step thickens the domain 
walls; the width of the domain wall near the step is about three times that of the prefect 
Pb-centered domain wall. Note that domain wall thickness cannot be defined uniquely very 
much like boundary layer thickness in fiuid mechanics. Here, domain wall thickness is by 
definition the region that is affected by the domain wall, i.e. those layers of atoms that are 
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distorted. One can use definitions like the 99%-thickness in fluid mechanics and deflne the 
domain wall thickness as the length of the region that has 99% of the far fleld rigid translation 
displacement. What is important here is that no matter what deflnition is chosen, domain 
wall "thickness" increases by the presence of steps. Note that due to the symmetry of the 
Pb-Pb step, atomic conflguration for negative values of s and d will have the same behavior. 
Also note that the y-components of the atoms on the section s = a in Figj3]^a) are not 
symmetric because of the way we deflne this section (see Figj2]). As the coordinates of cores 
and shells are close to each other, we only plot the results for cores. Also because other 
types of atoms display a similar behavior, we do not plot their coordinates here. 

We follow Meyer and Vanderbilt to calculate the polarization proflle for each step. 
The polarization of unit cell i is calculated as 



of the unit cell i from the ideal lattice site. Wj denotes the weight for atom j. For example, 
for a Ti-centered unit cell we have Wxi = 1, wq = 1/2, and wph = 1/8. We have plotted the 
polarization Pj of the rows of unit cells on the section v-v with d = a (see Figj2]) in a Pb-Pb 
step in FigJD^a), where Pj = Pj/lPtl, with IP^I denoting the norm of the polarization of the 
bulk. We obtain the bulk polarization of 80.1 fiC cm~^, which is close to the published values 
81.0 jj,C cm~'^ {22! and 81.2 jjC cm~'^ j^. We observe that near the step, the domain wall has 
a mixed Bloch-Neel character. Denoting the polarization components by P = {Px,Py,Pz), 
where Px, Py, and Pz are polarization components in (100), (010), and (OOl)-directions, 
respectively, we observe that the polarization vector rotates out of the (OOl)-plane with 
the Bloch angle Ob = tan~i [P^/Py] (see Figg]). The maximum rotation angle as for Pb- 
centered domain wall is ~ 7.0°. Also the polarization rotates in the (OOl)-plane with 
the Neel angle ajy = tan~^ (^Px/Py). The maximum value of ajy for Pb-centered wall is 



ajy — 9.9° (compare this with = 1.43° for the perfect domain wall [13|]). The maximum 
value of the polarization in the (100) and (OOl)-directions are about 13.9% and 12.2% of the 
bulk polarization, respectively. 

Finally, we calculate the energy of the Pb-Pb step, Eph^ph. Similar to the domain wall 
energy, we deflne the step energy to be the difference in energies of the unit cells inside the 




(7) 



where e is the electron charge, fic is the volume of the unit cell, Z* is the Born effective 
charge tensor of the cubic PbTiOs bulk, and u* denotes the displacement of the j-th atom 
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computational box that are located on the domain wall with the step and bulk energy of the 
same number of unit cells, divided by the total area of domain wall in the system. This way 
we obtain the Pb-Pb step energy to be 157 mJm~^, which is greater than the Pb-centered 
domain wall energy that is 132 mJm~^ ^. 

Ti-Ti step: Figsl2](c) and (d) depict the y-coordinates of Ti cores in a Ti-Ti step for 
different sections v-v and h-h. Again because of symmetry, we plot the results only for 
positive values of s and d and also similar to the Pb-Pb step, by definition of the section 
s = a, the y-components of the atoms on this section in FigJSl^c) are not symmetric. Similar 
to the Pb-Pb step, we observe that the Ti-Ti step is localized, i.e., atomic distortions are 
confined to a 9a x 18a box in the (OlO)-plane. Again we observe that the step thickens the 
Ti-centered domain wall; the thickness of the defective wall is about three times that of the 
perfect domain wall. 

As is shown in FigJD^b), polarization has a mixed Bloch-Neel character near the step. 
For the Ti-Ti step, the maximum value of the Bloch and Neel rotation angles are as — 9.5° 
and ttN — 8.1° (compare this with aN = 1-0° in the perfect domain wall [13|), respectively. 
The maximum value of the polarization in the (100) and (OOl)-directions are about 8.1% 
and 16.8% of the bulk polarization, respectively. The energy of the Ti-Ti step is Exi-n = 
172 mJm~'^, which is larger than the energy of the Pb-Pb step and the energy of the 
Ti-centered domain wall, which is 169 mJm~'^ jsl. This is consistent with the fact that 
Ti-centered 180° domain walls have a greater static energy than Pb-centered domain walls 

aa. 

Pb-Ti step: We have plotted the y-coordinates of Pb cores in a Pb-Ti step in FigJS] 
for different v-v and h-h sections. Note that because Pb-Ti steps are not symmetric, we 
have plotted the results for both positive and negative values of s and d. Also since other 
types of cores and shells have a similar behavior, we do not plot their coordinates here. We 
observe that similar to the other two steps, the Pb-Ti step causes local distortions that are 
confined to a 6a x 16a box in (OlO)-plane. The step broadens the domain wall; the defective 
domain wall thickness is twice that of the perfect domain wall. Note that the Pb-centered 
and Ti-centered domain walls for this step are half a lattice spacing apart and this may 
explain the weaker thickening effect of the Pb-Ti step. 

As Figj6] shows polarization distribution has a mixed Bloch-Neel character near the Pb-Ti 
step but the Bloch character is more dominant. The Polarization profile is plotted for the 



9 



row of unit cells on the section v-v with d — a, which is located in the Ti-centered part 
of the step. For the Pb-Ti step, the maximum value of the Bloch and Neel rotations are 
aB — 23.0° and ajy — 5.9°, respectively. The maximum value of the polarization in the (100) 
and (OOl)-directions are about 5.3% and 20.5% of the bulk polarization, respectively. The 
energy of the Pb-Ti step is Epi,_Ti = 165 mJm~'^. It is seen that Eph_pi, < Eph_Ti < Exi-n- 

V. CONCLUDING REMARKS 

In this work we obtained the atomic structure of three different types of steps on 180° 
domain walls in PbTiOs using the method of anharmonic lattice statics. We observe that 
these steps cause local atomic distortions that are confined to a box with dimensions of a few 
lattice spacings. All the three steps have a broadening effect on the domain wall thickness. 
Pb-Ti steps have a less broadening effect compared to the other two steps. 

We also observe that steps on 180° domain walls can cause the polarization profile to 
have a mixed Bloch-Neel character. The Bloch character is more dominant in Ti-Ti and 
Pb-Ti steps. Finally, we observe that the Pb-Pb step has a lower static energy than the 
other two steps. 
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Pb • Pb-centered domain wall 

Ti Ti-centered domain wall 

Origin * 

FIG. 2: A representative atomic layer for the initial configuration of the three steps: (a) Pb cores 
in Pb-Pb step, (b) Ti cores in Ti-Ti step, (c) Pb cores in Pb-Ti step, and (d) Ti cores in Pb-Ti 
step. Note that planes h-h and v-v are sections that are used for a better display of the variation 
of the distortion field in our numerical examples, s and d denote the distances of sections h-h 
and v-v from the reference planes and r^, respectively, and are parallel to (100) and 
(OOl)-planes, respectively, h-h and v-v sections in part (a) correspond to s = a and d = a and 
h-h and v-v sections in part (b) correspond to s = 2a and d = 2a. The shaded regions denote 
the computational box, which contains W x L unit cells and different colors show the regions with 
opposite polarization inside the computational box. The symbol in these figures denote the 
origin of the coordinate system in each type of steps. The blue and red filled and hollow circles 
denote the atoms whose displacements are used as the displacements of the atoms located outside 
of the computational box. 
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FIG. 3: The y-coordinates of atoms, (a) and (b) are Pb cores in a Pb-Pb step, (c) and (d) are 
Ti cores in a Ti-Ti step. Note that as it is shown in FiglJl s and d denote the distances from the 
reference planes. 
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FIG. 4: The polarization vectors P = {Px,Py,Pz) for the row of unit cells on the section v-v with 
d = a for: (a) Pb-Pb step and (b) Ti-Ti step. Close to the step, polarization rotates out of the 
(OOl)-plane with the Bloch angle as- The polarization also rotates inside the (OOl)-plane with the 
Neel angle a^- Note that the Bloch and Neel components of the polarization correspond to the 
components in (OOl)-direction [Pz) and (lOO)-direction {Px)^ respectively. 
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FIG. 5: The y-coordinates of Pb cores in a Pb-Ti step. Note that as it is shown in Figj2l s and d 
denote the distances from the reference planes. 
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FIG. 6: The polarization vectors P = {Px,Py,Pz) for the row of unit cells on the section v-v 
with d = a for a Pb-Ti step. Close to the step, polarization rotates out of the (OOl)-plane with 
the Bloch angle as- The polarization also rotates inside the (OOl)-plane with the Neel angle oat. 
Note that the Bloch and Neel components of the polarization correspond to the components in 
(OOl)-direction (P^) and (lOO)-direction (Pr), respectively. 
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